9 On Transformed Hamiltonians 



9.1 General 

If the four-component Dirac equation is approximated by some two-component 
form, certain savings can be obtained. Typically, the small components are 
eliminated and two separate equations are introduced for atomic j = I ± ^ 
states. For molecules, one typically divides the Hamiltonian into a quasirel- 
ativistic (spin-orbit averaged) part and the spin-orbit part, diagonalized at 
some stage. Special attention must be given to variational stability. 

One can follow these developments in Table 2.2 of RTAM I-III and the 
applications in Table 7.3 of RTAM III. For reviews, see Hess and Marian 
(2000a). 



9.2 The Foldy-Wouthuysen transformation 

Foldy and Wouthuysen (1950) presented a systematic procedure for de- 
coupling the large and small components of the Dirac equation, 

ho = ca- {p- eA) + I3mc^ + e$ (9.1) 
= f3mc^ + + 



(with $ the Coulomb field) where the "odd", 0, and "even", £, parts satisfy 

(9.2) 



= ca-{p-eA), {/3,0} = 0, 
£ = e$, [13, E] = 



Any operator H can be separated to 



0= ]^{H-(3H(5) 
£= '^-{H + (5H(5) 



(9.3) 



We shall annihilate the " odd" part up to some power of 1/ mc^ through the 
transformation 

= e^'^* (9.4) 



A good choice is 



S=^ = -i^^ (9.5) 
2mc^ 2mc 



where vr is represents the generalized momentum. 
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A third order FW Hamiltonian becomes 



h'n =mc^?e^ + — vr^ (9.6a) 
2m 

- —ps ■ B (9.6b) 
m 

— - — iri^p^ (mass-velocity) (9.6c) 
+ TT^^ ■ ^ (Darwin) (9.6d) 

:r^S_ -Exp (spin-orbit) (9.6e) 

2m-^c/ - 

eh 1 

„ o o iS- B)Tr^ (rel. corr. to Zeeman) (9.6f) 

m 2m^c^ 

e is the charge of the particele, h= 1, s = \a. [See Bjorken-Drell (1964), 
eq. (4.5), Moss (1973), eq. (9.75)]. 

This expression can be used in Ist-order perturbation theory like the 
Pauli one. It is not useful in variational calculations as both the mass- 
velocity term hm = ~^P'^ the Darwin term ho = ^'^K — —^-S{r) 

are strongly singular. One could just as well use them to first order, 0(c~^). 
Contracted basis sets could however allow variational use. 



9.3 The Cowan-Griffin equation 

Cowan and GrifRn (1976) or Wood and Boring (1978) as well as other 
authors use the 2nd-order equation of Dirac (1928a), putting all the norm 
in the large component. 



P'+-P+ i-mc + ^-^)Q = 
r c 

Q' --Q + {-mc+^^)P = {) 
r c 



(9.7) 



Q=^--^— — (P' + -P) (9.8) 
mc"^ + E — V r 



O' = - ip" +-P' - —p) + — (p' + 1)p (g g) 

jnc^ + E-V^ r r^"^' {mc^ + E - Vf^ r'^ ' 
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:{P" + ^P' - - - , " ^ (P' + -P) 



mc^ -\- E — V r -j-^ ' ^ j^c^ + E — V 

(mc^ + E — Vy ^ r ' c 



P"+-P' - -,P--P' - + -^^ [P' + -P) 

+ \{2mc^ + e-V){e-V)P = 



2 



1 



+ -^{e-Vy = (9.10) 



k(/« + 1) = /(l + l) (9.11) 

Quasirelativistic usage: put P^dr = l! (This is the only approx- 
imation!) 

Dirac (1928) [PRSL A 117 (1928) 610]: Exact 2°'i-order equation. 
For the non-relativistic case, 



r^oo: P" = -2m£P ^ P oc ^ = V-2m£ (9.13) 



r^O: P"='^^^P ^Par'+i, R o: ^ = (9.12) 



For the relativistic case, 



c 



9.4 Douglas-Kroll-Hess 

Instead of expanding in powers of 1/c one can rather expand in powers of 
the coupUng strength, Zach (Douglas and KroU 1974, Hess 1986, Jansen 
and Hess 1989). This leads to operators which are bound from below and 
can be used variationally. 

The DK formalism, which is an all electron method, is based on a series 
of unitary transformations Uo,Ui,. . . of which the lowest is the free-particle 
FW-trans formation defined by 

f Uo = A{1 + ^R) 

I , (9.15) 
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where we have 



A - (9.16) 



ca ■ p 
Ep + mc^ 

Ep = cVp^TmV (9.18) 



Applying Uq to ho = ca • P + — l)mc^ + V gives 

UohoUQ^ = pEp + S,i + Oi = Hi (9.19) 

with 



£i = A{V + RVR)A 
Oi = PA{RV - VR)A 

For chemical purposes one more transformation is needed. One uses 



(9.20) 



Ui = ^l + Wl + Wi (9.21) 

with W^i anti-hermitian, wl = —Wi. Performing the transformation through 
Ui and expanding the square root in powers of Wi, 

UiHiU^^ = f3Ep - [f3Ep, Wl] + £i + Oi 

+ ^pEpWf + ^WfpEp - WipEpWi 

+[W^i,Oi] + [H^i,£i] + ... (9.22) 

omitting higher order terms. The first-order odd term is eliminated by 
setting 

[pEp, Wl] = Oi (9.23) 
and solving for Wi. It's momentum space integral operator is 

Wi<^{p) = J (fp'Wi{p,p^)<^{^) (9.24) 

with a kernel 

W,{p,p') = AiR-R')A'^^^ (9.25) 

V{p,p') is the Fourier transform of the external potential. The primed quan- 
tities are expressed in terms of the variable p'. The final result is 

^decoupled ^pEp + Ei-P [WiEpWi + I [Wl Ep]] . (9.26) 
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The spin-orbit terms can be separated using 

{a- A){a- B) = A- B + ia- {Ax B) (9.27) 
see Hess et al. (1995). The final result becomes 



cl. IIUC. . , 

TT+ \ " \ " 2ry rZ-ia ^ 



' Ei + m Via Ei + m 
y\ -^(^ X p.) . {a, + 2a,)#^ (9.28) 



A series expansion in powers of c ^ becomes 
A 1 3p2 



Ei + mc^ 2mc^ 16m^c^ 



+ . . . (9.29) 



The DK Hamiltonian has now been implemented in a large number of 
programs, to mention a few: MAGIC, MOLCAS, PARAGAUSS and those by 
Hess himself. 

For calculating properties, the operator has to be transformed as well, if 
transformed Hamiltonians are used. This could be quite tedious. 

9.5 Zero Order Regular Approximation, ZORA 

The previous elimination procedures treated as a small parameter. 

Sufficiently near the nucleus this is not true. Here 

. . / V -E-\ ^/V- E\n 



n=0 



will only converge when |^^| < 1- By rewriting the expression as 
we obtain an expression in 



2mc2 - V 

which will converge. This leads to the Zero (or higher) Order Regular Ap- 
proximation, ZORA (E. van Lenthe et al. 1993-). 

The ZORA-level approximation was already given by Chang, Pelissier 
and Durand (1986) or by HeuUy et al. (1986). 

One problem with ZORA is the electric gauge invariance but this can 
be approximately remedied (E. van Lenthe, Baerends, Snijders 1994, van 
Wiillen 1998). 
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The scalar scaled ZORA equation is (Faas et al. 2000) 

\\p-Bp + Vc- K^^] = er'^'Ui (9.32) 



S 

with the scaling factor 

1 



S = l + {^i\p- ^B'pl^i) (9.33) 



c 



c2 



J and K are the usual Coulomb and exchange operators, respectively. Van 
Wiillen's idea was to set B = Bq = (? j (2c^ — Vb), Vb depending on the atoms 
only. 

ZORA has been implemented in packages as GAMESS-uk, adf and dirac. 
F. Wang et al. (2000) include spin-orbit effects at DFT level in 

The separated SO operator becomes 

hso = i^-[px ^^^J'_y^^ p\ (9.36) 

9.6 Direct Perturbation Theory, DPT 

The wave function and its metric will have to be an analytic function at the 
limit — > 0. 

The Direct Perturbation Theory (DPT) has been developed by Sewell 
(1949), Titchmarsh (1962), Rutkowski (1986-), Kutzelnigg (1989-). 
See also Franke (1994). 

Instead of the original Dirac equation 



one now writes 



V-E ca-p I I V'+ 
ca ■ p -2mc^ + V - EJ \V- 



V-E a-p i I 

a-p -2m + t{V - E)/c^J \cip. 



(9.37) 



(9.38) 
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with the parameter t = and t = 1 corresponding to the non-relativistic and 
relativistic cases respectively. The NR case corresponds to the Levy-Leblond 
(1967) equation. Both E and 




are expanded in powers of t (or c ^). A key paper on the subject is Kutzel- 
nigg (1996). 

9.7 Further examples 

9.7.1 RESC 

Hirao (see Nakajima et al. 1999, 2000) start from the inconvenience of 
having the energy in the 2nd-order equation 

+ • P) 2^,2 'V + E • ^)] V'L = £^^L (9.39) 

(V'lII + x^xIV'l) = 1, V's = X'4}i^ 
They suggest replacing 

E-V ^T= Vm2c4 + p2c2 _ mc^ (9.4O) 

which is independent of the eigenvalue E. The Hamiltonian can again be 
separated into spin-free and spin-dependent parts. 

Barysz (2000) has analyzed this family of RESC (Relativistic scheme 
for Eliminating Small Components) Hamiltonians. Although they can show 
signs of variational collapse already at moderate Z {k, 80), the result for 
valence electrons can be close to Douglas-KroU ones. 

9.7.2 AMFI 

The Atomic-Mean-Field-Integral method eliminates the spin-othcr-orbit terms 
in terms of an effective one-electron, one-centre potentials. (Hess et al. 1996, 
Marian and Wahlgren 1996, Schimmelpfennig et al. 1998ab, Tatchen and 
Marian 1999). 

An average over the two-electron terms is taken for the valence shell, 
and thus large two-electron integral files can be avoided. 
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